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1. INTRODUCTION 

Optimal reactive power problem has been key problem in power system, since it plays major role in 
secure & economic operation of the power system. Many conventional methods [1-8] have been applied for 
solving optimal reactive power problem. But many drawbacks have been found in the conventional methods 
and mainly difficulty in handling the inequality constraints. Last two decades many evolutionary 
algorithms [9-18] continuously applied to solve the problem. This paper proposes a new Sperm Motility 
(SM) algorithm for solving optimal reactive power problem. Sperm Motility (SM) algorithm is stimulated by 
fertilization process in human beings [19]. Wolf optimization (WO) algorithm has been formulated on the 
basis of basic deeds of wolf search for the prey. The formulation has been enhanced by utilizing the velocity 
& movement properties of particle swarm optimization algorithm.Both the projected algorithms have been 
tested in standard IEEE 57,118, 300 bus test systems. Simulated outcomes reveal about the reduction of real 
power loss & with variables are in the standard limits. Almost both algorithms solved the problem efficiently, 
yet wolf optimization has slight edge over the sperm motility algorithm in reducing the real power loss. 


2. PROBLEM FORMULATION 
The key objective of the reactive power problem is to minimize the system real power 
loss & given as, 
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Pioss= Zka , EKV +VÝ-2Vi Vi cos 04) ) 
Voltage deviation magnitudes (VD) is stated as follows, 

Minimize VD = YL, |V, — 1.0] (2) 
Load flow equality constraints: 
Poi- Ppi — V onb q | a rosee 0,i = 1,2 ...., nb (3) 
iXj=1 Vj [+B sin Bj 
Gj  sin®;j ; 
Qai — Qni - Visn.y, Fe cos 6 =0,i=1,2....,nb (4) 
Inequality constraints are: 
Vali" < Voi < Val", i € ng (5) 
ymin < Vu < Ve i e nl (6) 
ci S Qa s QE i Enc (7) 
GS Qai S QUI E ng (8) 
TEY < T; Te Te nit (9) 
smin < gmax j € nl (10) 


3. SPERM MOTILITY ALGORITHM 


Sperm Motility (SM) algorithm has been inspired by the fertilization procedure in human beings. 
Throughout the exploration sequence as the species chemo attractant, sperms are fascinated towards 
the ovum. Also Chemo attractant & concentration will induce the sperm when is moving closer to ovum. 
Utmost eminence sperm will be moved over & expressed as -Type A. With probability P, € [0,1] low quality 
sperms are discarded; specified as type B, C and D. Towards the ovum more than 220 million sperms swim 
capriciously with velocity v; at position x; , & by the Stokes equations motility can be described as, 


Re= (4 v.vv) + vp = uV’?v + fV.v=0xEN 
Simpler form of Stokes written as: 

Vp =uV*v + f 

V.v=0,y EN 


C singularity velocity solution as follows, 





vi(t) = (—) * (“H+ =) * Fly = (—) * sj, €) * Fly; i,j = 1,2,3... 


8ru h h3 8ru 


Due to a force Fl; concentrated at the point ¢, the flow will be as, 


h=y-¢ 
k? = h +h? +h 


It has been updated as, 


Yid = yE) + (Ë) «(Vier + vO) + a0 - 9 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 
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Chemo attractant is defined as follows, 
ca;(t) = cao (t) + ca (IIJ* — 0O (18) 
In current iteration J* is the outstanding solution exist. 


Sperm Motility(SM) algorithm for solving reactive power problem 
Commence 
Based on the problem objective function is defined. 
N sperm Population size is initialised 
Primary attentiveness co has been engendered for N sperm with reference to primary position yo and 
velocity vo are produced. 
Parameters of the motility are described; While (t< Maximum Generation) 
For i=1: N do 
Velocity vi is computed from By equation (14) using the data at t = ti; 
Position xi has been modernized for sperm i , by equation (17) 
Calculate the each sperm individual value, according to its position 
When new solution is available & superior then modernize the population. 
Value of cai from equation (18) has been calculated. 
When ca; < cai-1 , poorer sperm with help of (Pa), will be abandon 
Checking of the constraints with respective to objective function 
End for 
Existing outstanding population has been sorted out 
End 


4. WOLF OPTIMIZATION 

Wolf optimization mimics the communal management and hunt deeds of wolf in nature [20]. There 
are three fittest candidate solutions assumed as @,f and y to lead the population toward promising regions of 
the exploration space in each iteration of wolf optimization. g is named for the rest of wolves and it will 
assist a, and y to encircle, hunt, and attack prey, that is, to find Enriched solutions. In order to scientifically 
replicate the encompassing behavior of Red wolves, the following equations are proposed: 


= |Ë. Yp (©) — rO], 


G =| 
Y(t +1) =Y,(t) —H.G (19) 





Where t indicates the current iteration, H = 2b Ti —b, F= 2.T> , Yp the position vector of the 


prey, Yis the position vector of a wolf, b is linearly decreased from 2.0 to 0, and ryand ry are arbitrary 
vectors in [0, 1]. Hunting behavior of wolves are mathematically simulated by following equations, 


a= [F7 -7| 


Gp = |F Yg — ¥| (20) 
G, = |F3,¥, -Y| 
Yi = Yu — H; Ga 
Y, = Ý; — H, .Gg (21) 
Y; =Y, - Hz .G, 
Yi + 1) = 2 (22) 


3 


Position of wolf was updated by (19) & the following equation is used to discrete the position. 


1 Y,; > 0.50 


; (23) 
0 otherwise 


flagij = { 
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Where i, indicates the jth position of the ith wolf, flag; j is features of the wolf. 
To enhance the search velocity & position updating equations form particle swarm optimization has 
been incorporated in this approach. 


ving = Wp vi + cgy.Rm,.(mi — yi) + cgz. Rm3. (mÊ — yi) (24) 


Yi = Yet View (25) 


The current position of particle is yi & search velocity is vi . Global best-found position is. mË. In 
uniformly distributed interval (0, 1) Rm, & Rm, are arbitrary numbers. Where cg, and cg» are scaling 
parameters. is the particle inertia. The variable w is modernized as 


(tmax-0) 


tmax 


Or = (Wmax T, Omin): + Omin (26) 
Maximum and minimum of w+ is represented by Wmax and wWmin ; maximum number of iterations is 
given by tmax. Until termination conditions are met this process will be repeated. In this approach wolves 


categorized asa,fB and y determine the position of the prey. H= 2b .T; — b directs the exploration & 
exploitation process by reducing the value from 2 to 0.When |E] < 1 it converged towards the prey & If 


|H] > 1 diverged away. The first best Minimum loss and variables are accumulated as "a" position, score & 
as like second best, third best accumulated as "B" and " y" position & score. 

Commence 

Initialize the parameters 

Initialize b, H and F; beginning positions of wolves has been stimulated. 

i = 1: population size 

j=lin 

When (i, j) > 0.500 

@=1; 


End for 

Work out the maximum fitness of wolves as follows, 
Primary maximum fitness of the wolf is designated as “a” 
Second maximum fitness of the wolf is designated as “p” 
Third maximum fitness of the wolf is designated as “y” 
While k < maximum iteration 

For i= 1: population size 

Exact Location of the existing wolf has been revised periodically 
End for 

For i = 1: population size 

For i=1:n;If (i, j) > 0.496 

@e=l; 


End for 


Sporadically revise the values of b, H and F; 

At this stage Fitness of wolves has been calculated 

The assessment of wolves "a","B" and " y" has to be revised 
k=k+1; 

End while 

Re-examine the value of” æ “as the optimal characteristic division; 
End 


5. SIMULATED OUTCOME 
At first IEEE 57 bus system [21] is used as test system to validate the performance of the proposed 
algorithms. Total active and reactive power demands in the system are 1240.68 MW and 330.82 MVAR, 
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respectively. Generator data the system is given in Table 1. The optimum loss comparison is presented in 
Table 2. Number of iterations taken is 32 & time taken 11.24 sec. 


Table 1. Generator data 
Generator No Pgi minimum Pgi maximum Qgi minimum Qgi maximum 








1 25.00 50.00 0.00 0.00 

2 15.00 90.00 -17.00 50.00 
3 10.00 500.00 -10.00 60.00 
4 10.00 50.00 -8.00 25.00 
5 12.00 50.00 -140.00 200.00 
6 10.00 360.00 -3.00 9.00 

7 50.00 550.00 -50.00 155.00 





Table 2 comparison of losses 








Parameter Method Method Method Method Method Method Method SM WO 
CLPSO DE GSA OGSA SOA QODE CSA 
[22] [23] [23] [24] [22] [23] [25] 
PLOSS 24.5152 16.7857 23.4611 23.43 24.2654 15.8473 15.5149 12.1482 12.0064 
MW) 





Secondly IEEE 118 bus system [26] is used as test system to validate the performance of 
the proposed algorithms. Table 3 shows limit values and Table 4 show the comparison of results. 


Table 3. Limitation of reactive power sources 
Bus number 5 34 37 44 45 46 48 
Maximum value of QC 0.000 14.000 0.000 10.000 10.000 10.000 15.000 
Minimum value of QC -40.000 0.000 -25.000 0.000 0.000 0.000 0.000 
Bus number 74 79 82 83 105 107 110 
Maximum value of QC 12.000 20.000 20.000 10.000 20.000 6.000 6.000 
Minimum value of QC 0.000 0.000 0.000 0.000 0.000 0.000 0.000 











Table 4. Evaluation of results 








Active power loss — Minimum & Methodology - Methodology - Methodology SM WO 
Maximum values BBO ILSBBO/ ILSBBO/ 
[27] strategy 1 strategy 1 
[27] [27] 
Minimum value 128.770 126.980 124.780 126.540 = 125.172 
Maximum value 132.640 137.340 132.390 132.860 131.248 
Average value 130.210 130.370 129.220 129.120 128.864 





Finally IEEE 300 bus system [21] is used as test system to validate the performance of the proposed 
algorithms. Table 5 shows the comparison of real power loss obtained after optimization. 


Table 5 comparison of real power loss 





Parameter Method EGA [28] Method EEA [28] Method CSA [25] SM WO 





PLOSS (MW) 646.2998 650.6027 635.8942 630.1898 629.2824 





6. CONCLUSION 

In this paper both Sperm Motility (SM) algorithm & Wolf Optimization (WO) algorithm solved 
the problem successfully. Performance of both above said algorithms in solving the problem is outstanding. 
Both the projected algorithms have been tested in standard IEEE 57,118, 300 bus test systems. Simulated 
outcomes reveal about the reduction of real power loss & with variables are in the standard limits. Almost 
both algorithms solved the problem efficiently, yet wolf optimization has slight edge over the sperm motility 
algorithm in reducing the real power loss. 
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